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This  paper  presents  a  study  concerning  a  sensorless  vector  control  of  an 
induction  machine  fed  by  a  voltage  source  inverter.  The  aim  is  to  provide  a 
scheme  to  control  the  speed  and  the  rotor  flux  using  a  sensorless  integral 
backstepping  control  approach.  The  rotor  speed  estimation  is  done  by  an 
observer  using  the  model  reference  adaptive  system  (MRAS)  technique 
whereas  the  nonlinear  backstepping  observer  is  used  to  get  the  rotor  flux.  The 
main  objective  is  to  achieve  a  robust  control,  adaptive  and  efficient,  which 
will  allow  us  to  test  and  evaluate  the  performance  of  the  proposed  observer, 
combined  with  a  sensorless  control  of  the  induction  machine.  Tests  and 
validation  are  done  using  numerical  simulations  with 
MATLAB/SIMULINK-PSB  (Power  System  Block  set)  toolbox.  The  results 
show  good  performance  in  terms  of  robustness  regarding  machine  parameter 
variations  and  show  the  excellent  quality  of  the  control  law  associated  with 
the  observer,  despite  the  observability  problems  when  the  machine  operates 
at  low  speed. 
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1.  INTRODUCTION 

In  the  recent  decades,  the  variable  speed  electric  drives  developed  rapidly  and  aroused  significant 
interest  in  most  sectors  of  industrial  activity  implementing  the  driving  force,  due  to  technological  advances, 
but  also  methodological  [lj.The  direct  current  machine  has  long  been  the  best  source  of  electro  mechanical 
variable  speed  because  of  its  ease  of  control.  This  machine,  with  its  two  coils  to  orthogonal  fields,  presents  a 
natural  decoupling  between  stator  flux  and  electromagnetic  torque  [2].  However,  its  major  disadvantage  is 
the  existence  of  the  "collector-brush"  device  which  is  the  origin  of  his  employment  limitations.  This 
mechanical  device  is  poorly  tolerated  in  some  applications  for  its  high  price,  its  frequent  and  costly 
maintenance,  but  also  because  of  the  presence  of  electrical  arcs  prohibiting  use  in  explosive  atmospheres. 
Due  to  these  constraints,  research  was  carried  for  other  solutions,  exploiting  the  AC  machines  and  especially 
the  squirrel  cage  induction  motors  (IM)  [3]. The  IM  has  certain  advantages  related  to  its  simple  design, 
robustness,  low  cost  and  maintenance.  Also,  it  is  widely  used  in  many  fields  such  as  aerospace,  rail  and 
marine  traction,  transportation  systems  and  in  a  multitude  of  industrial  applications  (machine  tools,  driving 
machines,  pumping  systems,  ventilation  and  lifting  [4]. 

However,  the  use  of  the  induction  machine  for  variable  speed  drives  raises  a  couple  of  difficulties, 
which  leads  to  a  complex  control  compared  to  that  of  another  motor.  Indeed,  this  type  of  motor  is 
characterized  by  a  multivariable  and  nonlinear  mathematical  model,  with  a  strong  coupling  between  the 
control  variables  of  the  magnetic  flux  and  the  electromagnetic  torque.  It  is  not  possible  to  control  a  motor 
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output  variable  (torque,  speed  or  position)  by  varying  a  single  input  variable  (current  or  voltage).  Moreover, 
some  of  its  state  variables  are  not  accessible  to  measurement  and  the  motor  parameters  are  subject  to 
variations  in  time  [5],  [6]  .The  implementation  of  these  control  methods  were  only  possible  through  the 
development  of  digital  signal  controllers,  microcontrollers,  static  converters.  Also,  several  techniques  have 
been  proposed  to  obtain  for  the  induction  machine,  the  high-performance  control  devices,  and  experience  its 
qualities  and  benefits  in  the  field  of  variable  speed  [7]. 

Among  the  controls,  we  distinguish  the  vector  control  based  on  dynamic  models,  using  the  Park 
transformation  to  solve  the  torque-flux  decoupling  issue.  The  combination  of  this  technique  with  a  modern 
automatic  approaches  based  on  nonlinear  control  algorithms  such  as  the  Backstepping  control  contribute  to 
solve  the  robustness  problem  (the  change  of  load  torque  and  machine  parameters).  However,  it  requires  the 
installation  of  an  incremental  encoder  to  measure  the  rotor  speed  [8-9].  Using  the  encoder  leads  to  additional 
cost  which  can  be  higher  than  the  machine  one,  for  low  power  application.  In  addition,  there  must  be  an 
additional  space  for  installing  the  encoder.  The  system  reliability  decreases  because  of  this  fragile  device. 
Because  of  all  these  inconvenient,  one  tried  to  exclude  the  use  the  incremental  encoder  and  research  on  the 
sensorless  control  of  induction  machines  began.  Several  strategies  were  proposed  in  the  literature  to  achieve 
this  objective  [10-11].  Much  of  the  proposed  methods  are  based  on  observers  and  estimators  in  open  loop. 
They  largely  depend  on  the  model  of  the  induction  machine.  For  this  reason,  these  techniques  fail  to  replace 
the  incremental  encoder  in  the  field  of  low  speeds.  Further  research  based  on  the  combination  of  modern 
techniques  and  theories  including  adaptive  laws.  These  techniques  promise  to  give  better  results,  even  in  the 
low-speed  including  stopped. 

In  this  paper,  we  chose  to  use  a  nonlinear  integral  backstepping  control  design  scheme,  for  speed 
control  and  we  propose  a  new  combination  for  the  realization  of  robust  observers,  by  using  the  model 
reference  adaptive  system  (MRAS)  to  estimate  the  rotor  speed  and  the  Backstepping  observer  to  estimate  the 
rotor  flux. 


2.  INDUCTION  MOTOR  MODEL  AND  VECTOR  CONTROL 
2.1  Model  of  IM  in  (d-q)  Reference  Frame 

The  mathematical  model  of  the  induction  machine  in  the  (d-q)  synchronously  rotating  frame  is  given 
by  the  following  nonlinear  equations  [12-13]: 


With 
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2.2  Direct  rotor  flux  orientation  (DRFO) 

The  aim  of  vector  control  is  to  let  the  induction  machine  effectively  respond  to  reference  variations 
(rotor  position,  torque,  speed)  in  a  wide  range  of  operation  and  for  applications  requiring  high  dynamic 
performance.  It  is  based  on  the  instantaneous  torque  control  and  a  suitable  choice  of  reference  (d,  q)  of  Park, 
so  that  the  d-axis  is  aligned  with  the  rotor  flux  vector,  to  obtain  the  decoupling  between  the  electromagnetic 
torque  and  rotor  flux.  This  alignment  imposes  the  following  condition:  [14-15]: 

(pr  —  (prc[  and  tppq  —  0  (2) 


This  choice  allows  to  control  the  torque  via  the  current  iSq  and  the  rotor  flux  through  the  current 
^sd 
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3.1  Principle  of  the  integral  backstepping  control 

The  Backstepping  technique  is  a  relatively  new  control  method  of  nonlinear  systems.  The  principle 
of  this  approach  is  to  stabilize  at  the  beginning  the  first  sub  system  by  a  stabilizing  function  known  via  a 
chosen  Lyapunov  function,  then  to  add  to  its  error  an  integrator.  The  same  procedure  for  the  successive 
subsystems  to  finally  reach  a  global  Lyapunov  function  giving  the  overall  control  law  that  stabilizes  the 
system  [16-17].  The  combination  of  technical  Backstepping  with  vector  control  of  the  induction  machine 
gives  the  better  performance  of  the  control  law  [18-19]. 

3.2  The  integral  Backstepping  control  for  IM  and  rotor  flux  estimator 

The  objective  to  design  an  integral  backstepping  control,  allowing  to  ensure  the  tracking  of  rotor 
speed  and  rotor  flux  for  the  machine  [20-21],  we  define  the  errors  and  respectively  representing  the 

error  between  the  actual  speed  (Dr  and  the  reference  speed  CO*  and  the  error  between  the  flux  module  cp^. 
and  its  reference  cp^  . 

a)  First  step  "speed  and  flux  loop'' 


e(0,  =cor-cor+z(i 


e<p.  =<Prd  -<Prd  +Z, 


<Prf 


(3) 


With: 


Then  the  dynamical  error  equations  are  given  by: 


^cprd  ^  ^Prd  ^cprd 

With: 

=k;„r(M;-coi.) 

Then  the  error  dynamical  equations  becomes: 

^=®LMV^)+a7-wr+aA+\ 

V  =(Prd  “Msd  +a5,(Prd 


We  require  that  er,  =— k  e,  and  em  =— em  ,  where  k  and  k  are  positives  constants. 

1  CO,.  03,.  Cprd  (prd  Cprd  ®r  <Pdr  r 

We  can  use  the  Lyapunov  function: 


V 


(6) 
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Its  derivative  is: 


v=e  e  +e  e 

»  >rd^(p, 


v=eo,(H-a6.(isq.(pld)+a7.0r+a8.Tl+z()t)+e(pj(j);-a4isd+a5.(pld+ziP( 


(7) 


The  aim  of  tracking  is  fulfilled  ( V<0)  choosing  the  virtual  control  that  represents  the  stabilization  function  as 
follows: 


Cj  - [ke®,  + ' ®r  +  a7'®r  +  VT1  +  \  ] 

(Pid-a6 

a4 


The  derivative  of  the  Lyapunov  function  becomes: 


v  =  -kmrerar  -k(ftje9d  <0 


So,  the  errors  in  (3)  converge  to  zero. 

The  state  cprd  is  unknown,  so  we  need  to  use  his  estimation  in  (8),  so  the  new  expression: 


(8) 


+a7.cor  +a8.Tj  +zj 

a4 

b)  Second  step  "currents  loop" 

Once  virtuals  inputs  are  calculated,  we  define  the  error  on  the  currents  as  follows: 


eisq  Xsq  Xsq  +  ZiSq 


eisH  ~~  1sd  1sd  +  Zisd 


With: 


z- 


Z- 


Then,  the  error  dynamical  equations  are  expressed  by: 


eL  =  t — [k^e  +©;  +a7.w  +a8.T, +z  ]-i  +z 

sq  n\  n  r  r  r  n  sc 


Then,  the  error  equations  (5)  can  be  expressed  as: 

e»r  +a6-9rd-elsq  -a6.cprd-\ 

_  ~~kdr  e,Pcl,  +a4-eisd  -a5-9ri  a4'Zisd 


(10) 
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Where  the  flux  estimation  error  is: 

9rd=9rd-<Prd 

As  well,  the  dynamical  equations  for  the  error  ei  ,  ei  d  can  be  expressed  as: 


eisq  lsq  1sq  +  Zis 


<4  =  Vl,  ~K  bVsq  +  4 


Asq  Asq 


And 


^isd  1sd  +  Zisd 


Where: 


¥i  —— - [k^e  +6),  +a7.cor  +a8.T|  +z  ]  +  T - [k^e  +4  +a7.d)r  +a8.T,  +z  ] 

<Pri-a6 

Vi,=^K^+(Prd+a  A  +  4J 

a4 

\  =-®Sis6  +ai-iSq  -a3.(pri-®r  +a2-(PKI 
K  =aikd  +  ®s4,  +a2-tPrd  +a3®r¥rq 


Now  we  extend  the  Lyapunov  function  in  (6)  to  include  the  new  states  variables  ei 
parameter  estimation  errors  cprd  as: 


v=-(ei+e®  +ei2  +ei2  +-C) 

2V  COr  'Prd  'sq  'sd  y  '  rcl ' 


Where  y  is  a  positive  constant  known  as  the  adaptive  gain.  We  use  the  energy  function 

both  the  control  and  estimation  algorithms. 

The  derivative  of  the  above  equation  is: 


. L  , 

v  =e  e  +e  e  +e-  e  +e-  e  +-q),(d, 

^  ®r  9rd  <ftd  'sq  %  4d  'sd  y  ™  Yrd 


With: 
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,  e-  and  the 
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V ,  to  obtain 
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By  replacing  systems  (11),  (12)  and  (13)  in  (15)  equation  is  obtained: 


v  =-k  e2  -k  e2  -k  e2  -k  e2  +-(p2d+e.  (a^cLe., +k.  e  +M  -k  -bv+z. 

®r  ®r  9rd  ‘Prd  ’sq  ‘sq  lsd  'sd  y  ™  xsq  V  0  Tra  ®r  Xsq  Xsq  T  Xsq  Xsq  SCl  Xsq 


+%^+K%+%-K~h^+i0 


+9i-a5-efc  -r^e.  +-(%  -a^  +a5.q>J 


U5’>d  T  ^  v 

4n5  1 


For  the  condition  that  v  <  0 ,  the  control  voltages  Vsd  and  Vsq  as  follows: 

vsd  =  gL  (k-  e-  +w  -bvsd+z.  ) 

SCx  s  xsd  xsd  '  xsd  xsd  SCx  xsd  7 

v  =oLs(ki  e,  +v|/j  -X,  -bv  +Zj  ) 

^4  5  Asq  Asq  1  Asq  Asq  Asq 

A 

And  the  adaptation  law  of  rotor  flux  given  by: 

9rd  =-Y(-a5-ecprd  -rHd)  +  a4-isd  -as4d 

Mn5 

4.  ROTOR  SPEED  MRAS  OBSERVER 
4.1  Principle  of  MRAS  observer 

Currently  there  are  many  techniques  for  achieving  speed  estimator,  but  the  most  used  is  certainly  the 
approach  using  an  Model  Reference  Adaptive  System  (MRAS),  this  for  its  easier  implementation  and  its 
good  performance  in  terms  of  precision  and  stability  [22] . 

The  technique  of  MRAS  based  on  the  rotor  flux  for  estimating  the  speed  of  the  induction  machine, 
originally  developed  by  C.  Schauder  [23],  is  certainly  the  most  popular  strategy.  The  principle  of  the  speed 
estimation  by  a  MRAS  is  based  on  the  comparison  of  rotor  flux  outputs  of  two  estimators  using  different 
formulations.  The  first  estimator  is  based  on  the  stator  equations  defining  a  voltage  model,  called  the 
reference  model  and  the  second,  named  adjustable  or  adaptive  model  is  described  by  the  rotor  equations,  it 
defines  a  current  model. 

The  estimated  components  of  the  rotor  flux,  for  both  models,  are  based  only  on  the  measurement  of 
stator  quantities.  We  usere  the  reference  voltages  and  the  measured  currents  usually  expressed  in  a  stationary 
reference  frame  related  to  the  stator  (a-p)  [5].  The  structure  of  the  adaptive  observer  is  shown  in  Figure  1. 


(16) 


(17) 


Figure  1.  Block  diagram  of  MRAS  observer 
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4.2  Estimation  of  rotor  speed  by  MRAS  observer 

The  equations  of  the  reference  model  are  given  by: 

(Pra=7^(Vsa-RSisa-CTLSlsa) 

Lm 

< 

4>rp=7^(vsp-Rsisp-oLsisp) 

Lm 

The  equations  of  the  adjustable  model  are  given  by: 

9ra=“9ra-^9rp+Y'isa 
1  «  -  L  . 

(prp=-^r9rp+®r9ra+YLlsp 


(18) 


(19) 


With: 


We  notice  that  the  rotor  speed  is  only  present  in  equations  (19)  and  not  in  equations  (18).  The  adaptation 
algorithm  is  chosen  to  make  the  adjustable  model  converge  toward  the  reference  model  by  minimizing  the 
error.  For  this,  the  algorithm  parameters  are  defined  according  to  the  criterion  of  hyperstability  of  Popov. 
[23] 

The  error  between  the  states  of  the  two  models  can  be  expressed  in  matrix  form  by: 


<Pra  -9ra 

And  from  equations  (18)  and  (19)  its  dynamics  are  given  by: 

£=As-W 


(20) 


(21) 


With: 


1 


8  = 


,W=(cor-cor). 


Schauder  [24]  offers  an  adaptation  law  that  meets  the  criterion  of  Popov  and  given  by  equation: 


©r  =X2(s)+ j0Xi(s)dx 


(22) 


The  criterion  of  Popov  requires  satisfaction  of  the  following  inequality  [23]: 

JVwdt>-|u2 


(23) 


IJPEDS  Vol.  8,  No.  4,  December  2017:  1650-1662 


IJPEDS 


ISSN:  2088-8694 


□  1657 


With  |U  is  a  positive  constant. 

Replacing  8  and  W  by  their  values,  Equation  (23)  becomes: 


[{[yPrp  -£p9ra]  ®r  " %2 (S)  —  Jq Xl (S)dx 


We  can  solve  the  equation  (24)  using  the  following  relationship: 

{‘k(|f(t))f(t)dt>-3kf(0)2 

The  use  of  the  above  expression  provides  the  following  functions: 

Xi  =k1((£p9ra  -^Ap) =K(%9m  -«) 

h =kP(Sp9m  -s„q>rp) =kp(cprp9ra  -qyftp) 

Replacing  the  equations  (26)  in  equation  (22)  there  was  obtained: 
®r  =kp(9rp9,a  -9ra9rp)+ki  j(9rP9,a  "Mp# 


(24) 


(25) 


(26) 


(27) 


kp  andkj  are  the  PI  controller  gains  of  the  adaptation  mechanism. 
The  flux  position  and  the  Park  angle  is  determined  by: 


0S  =  arctg 


/  -  A 
9rp 

V^Pray 


(28) 


The  integral  Backstepping  control  with  the  MRAS-Backstepping  observer  combination  to  estimate  the  rotor 
flux  and  the  rotor  speed  for  induction  machine  can  be  presented  by  the  block  diagram  of  Figure  2. 


Figure  2.  Block  diagram  of  Sensorless  integral  Backstepping  control  scheme  associated  with 

MRAS-Backstepping  observer 


5.  SIMULATION  RESULTS  AND  INTERPRETATIONS 
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5.1  Simulation  results 

To  evaluate  the  performance  of  the  control  and  the  two  observers  proposed  in  this  study,  two 
Benchmarks  are  proposed  by  [3].  The  first  Benchmark  presents  rapid  transitions  and  operating  areas  at  large 
and  low  speeds,  the  second  Benchmark  consists  of  low  and  very  low  speed  operation  and  inversion. 

The  objective  is  to  impose  reference  trajectories  taking  into  account  the  problems  of  observability  of 
induction  machine  at  very  low  speed.  Recall  that  the  sufficient  condition  of  observability  loss  corresponds  to 
having  simultaneously  constant  mechanical  speed  and  zero  stator  pulsation  [25]. 

The  results  presented  below  were  obtained  by  simulation  tests  conducted  for  sensorless  integral 
Backstepping  control  scheme  associated  with  the  MRAS  and  Backstepping  observer,  on  the  profiles  imposed 
by  Benchmarks  1  and  2. 

a)  Benchmarks  "1" 
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Figure  3.  Performance  in  Benchmarks  1  of  the  proposed  observer  during  starting  operation  with  load  change 


b)  Benchmarks  "2” 
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Figure  4. Performance  in  Benchmarks  2  of  the  proposed  observer  at  low  speed  region 

5.2  Interpretations 

Benchmark  1  defines  the  reference  trajectories  for  speed,  rotor  flux  and  load  torque,  as  shown  in  the 
figure  3.  At  the  beginning,  reference  speed  is  imposed  zero,  to  enable  the  flux  to  reach  its  nominal  value  of 
lWb  which  is  then  kept  constant.  At  t=0,2  s,  the  rotor  speed  is  increased  to  100  rad/s,  according  to  an 
acceleration  ramp,  and  remains  constant  up  to  1,5  s.  Then  it  is  reduced,  reaching  a  value  of  zero  at  1,7  s,  the 
machine  is  then  maintained  at  standstill  until  2,5  s.  This  first  phase  is  used  to  test  the  behavior  of  the 
sensorless  control  for  speed  variations  in  quick  transitions,  and  for  a  first  critical  area  of  operation,  between 
1,7s  and  2,5s,  at  zero  speed  and  without  mechanical  load.  The  application  and  removal  of  a  torque  load  at  the 
moments  0,8s  and  1,2s  will  also  assess  the  impact  of  this  type  of  disturbance.  Then,  between  2,5  s  and  4,8  s,  a 
quasi-symmetric  velocity  profile  is  imposed  in  the  opposite  rotational  direction,  defining  a  second  area  of 
critical  operating  at  a  rate  of  -3,25  rad/s  between  4  s  and  4,8  s.  Finally,  the  machine  speed  is  again  increased 
to  100  rad/s,  to  be  there  maintained  up  to  6  s.  During  the  second  critical  phase,  a  torque  load  equal  to  50%  of 
the  nominal  value  is  applied  at  time  3,25  s  and  canceled  at  5,5s,  with  the  aim  of  analyzing  the  behavior  of  the 
training  in  an  unobservable  area,  but  this  time  in  the  presence  of  a  mechanical  load.  The  objective  of 
Benchmark  2  is  to  evaluate  the  performance  strategies  of  sensorless  control  at  low  and  very  low  speed  sensor 
in  conditions  of  unobservable  area  of  the  machine  and  to  highlight  the  influence  of  the  load  torque  on  the 
behavior  of  the  control. 

In  figure  3,  we  notice  that  the  estimated  speed  perfectly  follows  the  machine  speed  with  good 
behavior  in  terms  of  tracking  and  disturbance  rejection.  The  tracking  error  is  minimal  never  exceed  1  rad/s 
during  transient.  The  estimated  speed  converges  quickly  and  accurately  to  the  actual  machine  speed  and  the 
estimation  error  is  also  small,  but  presents  some  oscillations  between  1,25  rad/s  to  1,40  rad/s,  during  motor 
start.  When  going  from  100  to  0  rad/sec,  then  -100  to  0  rad/s,  this  error  is  0,3  and  0,4  rad/s,  while  in  the  two 
critical  areas  the  maximum  values  are  0,050  and  -0,040  rad/s  for  steady  state.  The  decoupling  is  adequately 
maintained.  However,  we  observe  small  oscillations  during  transients  and  very  low  orientation  errors  in  the 
two  critical  areas. 

Figure  4  shows  clearly  that  the  estimated  rotor  speed  follows  with  excellent  precision  the  actual 
machine  speed.  The  estimation  error  is  very  small  at  low  and  very  low  speed  steady.  When  reversing  the 
direction  of  rotation  from  -8  rad/s  to  8  rad/s,  the  estimation  shows  some  oscillations  with  a  maximum  error  of 
0,4  rad/s,  which  is  acceptable. 


6.  CONCLUSION 

In  this  paper,  a  new  sensorless  drive  system  approach  is  presented,  coupled  with  a  direct  vector 
control  for  induction  machine.  The  proposed  method  consists  on  the  use  of  integral  backstepping  control  with 
a  nonlinear  observer  of  rotor  flux  and  adaptive  observer  for  estimating  rotor  speed. 

The  results  show  excellent  estimation  under  different  operating  conditions,  especially  at  low  speeds. 
In  addition,  it  shows  robustness  against  external  disturbances  and  changes  in  machine  parameters.  This  study 
shows  that  it  is  possible  to  realize  an  association  of  two  observers  for  sensorless  control  based  on  a  nonlinear 
approach  with  a  high  level  of  performance. 


7.  NOMENCLATURE 

Tem,TL :  Electromagnetic  and  load  torque. 
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isd,iSq :  Stator  currents. 

Cprd,cpK1 :  Rotor  flux. 

C0m :  Electrical  angular, 
a:  Total  leakage  factor. 

Rs  and  Rr :  Stator  and  rotor  resistances. 

Ls  and  Lr :  Stator  and  rotor  inductances. 

Lm :  Mutual  inductance, 
p :  Number  of  pole  pairs. 

J  :  Moment  of  inertia. 
fc  :  Friction  coefficient. 

8.  APPENDIX 

Table  1:  Nominal  parameters  of  the  induction  machine  used: 


Rated  power 

1.5kW 

Rated  voltage 

220V 

Rated  speed 

1428rpm 

Nominal  frequency 

50Hz 

Rated  stator  current 

3.64  A 

Stator  resistance 

4.85Q 

Rotor  resistance 

3.805Q 

Cyclic  stator  inductance 

0.274H 

Cyclic  rotor  inductance 

0.274H 

Mutual  inductance 

0.258H 

Number  of  pole  pairs 

2 

Moment  of  inertia 

0.031kg/m2 

Friction  coefficient 

0.001 14Nm.s/rd 
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